Static point-to-set correlations in glass-forming liquids 
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We analyze static point-to-set correlations in glass-forming liquids. The generic idea is to freeze the 
position of a set of particles in an equilibrium configuration and to perform sampling in the presence 
of this additional constraint. Qualitatively different geometries for the confining set of particles 
are considered and a detailed comparison of resulting static and dynamic correlation functions is 
performed. Our results reveal the existence of static spatial correlations not detected by conventional 
two-body correlators, which appear to be decoupled from, and shorter-ranged than, dynamical length 
scales characterizing dynamic heterogeneity. We find that the dynamics slows down dramatically 
under confinement, which suggests new ways to investigate the glass transition. Our results indicate 
that the geometry in which particles are randomly pinned is the best candidate to study static 
correlations. 
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The collective nature of the dynamics of supercooled 
liquids approaching the glass transition is well estab- 
lished PJ. A recent key advance was the study of multi- 
point dynamical correlation functions instead of tradi- 
tional two-point correlation functions such as the inter- 
mediate scattering function. A central outcome is the 
determination of a dynamical length scale increasing by 
a factor of about 5-10 when glass-formers are cooled 
from normal liquid conditions to temperatures around 
the glass temperature T g . This raises the question of 
the underlying existence of nontrivial structural corre- 
lations, which also grow when the glass transition is 
approached, and how they relate to dynamic ones 
However, just as two-time dynamic correlation functions 
do not detect directly dynamic heterogeneity, two-point 
density correlation functions (pair correlation functions) 
seem unable to capture the relevant structural correla- 
tions, as these functions only show a mild temperature 
dependence. There are proposals that three-body ori- 
cntational order parameters might give insight for some 
specific glass-formers 0, H[ , or that local geometric struc- 
tures might be significant 0, Q , but generic methods to 
detect static order still need to be devised. 

Recently, the idea emerged that some form of 'amor- 
phous order' should develop in viscous liquids, which 
could be detected through 'point-to-set' correlation func- 
tions 04H1- Point-to-set (PTS) correlations probe static 
multi-point correlations, since they are determined by 
fixing the position of a 'set' of k particles and measuring 
the probability to find a (k + l)th particle at position 
rfe+i. It can be hoped that if the geometry of the frozen 
set is well chosen, these multi-point functions yield spa- 
tial information without measuring how the correlation 
function depends on all its k + 1 arguments. For a spher- 
ical cavity of radius d, for instance, one expects to detect 
a change of physical behavior when d interferes with the 
relevant structural length scale 0. Although first moti- 
vated in the context of the random first order transition 
(RFOT) theory [l2j], the set-up is actually more general 
and does not rely on any hypothesis regarding the mi- 



croscopic nature of the measured correlations. The gen- 
erality of the approach thus strongly suggests that it is 
important to explore in detail these PTS correlations in 
different geometries for the set as well as for various mod- 
els of glassy systems. Although reminiscent of studies of 
glass- formers in confined geometries (l3j . we emphasize 
that PTS correlations probe genuine bulk correlations, 
with no contribution from an external substrate 0, EiJ ■ 
While difficult to implement for molecular systems, in- 
vestigations along the lines suggested in the present work 
could be performed in colloidal materials where it is pos- 
sible to freeze the position of a selected set of particles 
using for instance optical tweezers. 

In this work, we show that PTS correlations can be 
detected using a broad variety of qualitatively distinct 
geometries, see Fig. [TJ which all reveal information on 
static correlations not included in conventional pair cor- 
relations. As a first step, we present the results of com- 
puter simulations of a simple glass-former in which we 
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FIG. 1: Four qualitatively different geometries to investigate 
point-to-set correlations in glass-forming liquids. 
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investigate how the geometry of the pinned particles in- 
fluences the measured correlations at a single state point. 
Our results suggest that the cavity geometry considered 
in earlier studies might not be the optimal choice to study 
static correlations. Another surprising outcome of our 
analysis is that static lengthscales appear to be decoupled 
from, and smaller than, dynamical lengthscales charac- 
terizing dynamic heterogeneity in the mildly supercooled 
regime typically studied in computer simulations. 

For each of the geometries shown in Fig. Q] we first 
equilibrate a three dimensional bulk system composed of 
N particles. At t = wc permanently pin k particles 
(filled circles in Fig. [1]) whereas the remaining N ~k par- 
ticles (open circles) move as before. Since the pinned 
particles were chosen from the equilibrated fluid, the 
thermodynamics of the free particles is strictly unper- 
turbed, provided an average is performed over both ther- 
mal fluctuations and different realizations of the pinning 
disorder 0, [Hj]. For the wall shown in Fig. [TJi, the 
average density profile is for instance strictly constant, 
(p(z)} = po, with no layering. We use constant temper- 
ature molecular dynamics simulations to study a 50:50 
binary mixture of harmonic spheres with size ratio f .4, 
all particles having the same mass m [l6[. The unit of 
length is given by <r, the diameter of the small parti- 
cles, the unit of time by yjma 2 / 'e, where e is the in- 
teraction strength, and temperature in units of 10 e, 
setting the Boltzmann constant ks = 1.0. We work at 
fixed density p = 0.675. For these parameters fl7j . 
slow dynamics sets in when T < 10, a fit to a mode- 
coupling divergence yields T c « 5.2, and the harmonic 
spheres behave as quasi-hard spheres, as discussed in de- 
tail in Refs. [16| . Thus our system is a canonical model 
for studies of the glass transition phenomenon. As an- 
nounced, we present a comparative study of the various 
confinement shown in Fig. [T] for a single, moderately low 
temperature, T = 8 > T c , a temperature at which the in- 
termediate scattering function of the bulk system already 
shows two-step relaxation. 

In all cases we wish to answer the following question: 
How does the presence of a pinned set of particles af- 
fect the structure and dynamics of the remaining free 
particles? To quantify these effects we define two over- 
lap functions, akin to the collective and self intermediate 
scattering functions. The collective overlap reads 

Q(t) = ^(n i (tK(0))/^(n i (0)), (1) 

i i 

where the sum runs over the cells (of volume v rs 0.53 3 , 
comparable to the particle volume) of a cubic grid, and 
rii(t) £ {0, 1} is the occupation number of cell i at time 
t. We set rii = if cell i contains a pinned particle. The 
overlap Q(t) is close to 1 if configurations at time and t 
are similar, but Q(t) is unaffected by particle exchanges. 
The long-time limit of the overlap, Qoo = Q(t — >• oo), 
provides direct information on static correlations. We 



also define the single particle function 

QselfW =£K(*K(0)>/£>?(0)>, (2) 

i i 

where nf (t) = 1 if the same particle occupies the cell % at 
times and t, and nf(t) = otherwise. 

We now describe the various geometries of Fig. [1] (a) 
Particles outside a cubic cavity of linear size 2d are frozen. 
The overlaps are measured at the center of the cavity, 
using the 4 3 central cells to improve the statistics; 8 in- 
dependent realizations are studied for each d. A similar 
(spherical) geometry has been studied in [l(| Hi} - (b) 
Particles outside the range < z < 2A are frozen, such 
that the free particles are 'sandwiched' between two in- 
finite walls separated by a distance 2A. The overlap is 
averaged over cells located in the plane parallel to the 
walls in the middle of the sandwich; 10 independent re- 
alizations are studied for each A, with L x — L y = 13.68. 
(c) A finite fraction, c, of particles is randomly se- 
lected in the fluid, with a typical distance between them 
21 = c -1 / 3 [lg|. The overlap is averaged over all cells. 
By using a large isotropic system, L — 16.3, only few re- 
alizations (typically 2-3) are needed to yield statistically 
accurate results, (d) Particles in the semi-infinite space 
z < are frozen. The overlap is averaged over cells be- 
longing to planes parallel to the wall at distance z from 
it. This geometry has been studied in more detail (in- 
cluding different temperatures) in [17[. For (b) and (d) 
we also included a hard wall at the boundary between the 
wall(s) and the fluid to prevent free particles to penetrate 
the walls |14j |. To perform quantitative comparisons, we 
define a 'confining length' £, respectively as £ = d, A, 
I, or z: a smaller £ means stronger confinement. Physi- 
cally, £ represents for each geometry the typical distance 
between the point where the overlap is measured to the 
pinned set of frozen particles. Note finally that geometry 
(a) is peculiar since the number of confined particles is 
always finite, while it scales with system size and diverges 
in the thermodynamic limit in cases (b-d). 

In Fig. [5] we gather our results for the four geometries, 
both overlaps (JTJ and © , and various degrees of confine- 
ment. In all geometries the time correlation functions 
have a similar qualitative behavior. When £ — > oo, bulk 
behaviour of the overlap Q(t) is recovered, with a two- 
step decay, and at long times a relaxation to the random 
value, Q r and ~ 0.110595 ss pqv. When confinement in- 
creases, the time dependence of Q(t) slows down, while 
the long-time limit increases, Qoc(0 > Qrand- In prac- 
tive we extract Qoo(£) by fitting the long-time decay of 
Q{t) to stretched exponential form. The quality of the 
fit is very good, as examplified in Fig. [2}d. In order to 
see that such a fit does indeed allow to obtain Qoo(C) 
with high precision, we show in Fig. [3] the time depen- 
dence of Q{t) — Qrand on a logarithmic scale as well as the 
corresponding fits. Finally we note that the evolution of 
Qself (t) is similar, showing a dramatic slowing down with 
increasing confinement, but its long-time limit is always 
zero. 



3 




t 



FIG. 2: Time dependence of the overlaps Q(t) (full lines) and 
Qseif(t) (dashed lines) for T = 8.0. The panels correspond to 
the four geometries of Fig. [1] In (b) we have also included the 
results of fits with a stretched exponential to the Q(t) data 
for A = 2.25 and A = 3.25 (dashed-dotted lines). 



A comparison between both functions shows that 
Q(t) roughly reaches Qao when Q se if (i) approaches zero, 
i.e. essentially when all particles have escaped the posi- 
tion they occupy at t — 0. Thus, when Qoo > Q ra nd, non- 
random local density fluctuations persist even though 
particles diffuse and explore the available space. By mon- 
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FIG. 3: Time dependence of the overlap Q(t) — Q ran d (sym- 
bols) and the stretched exponential fits (lines) for the sand- 
wich geometry and different values of A. 



itoring the ^-dependence of Qoo, we can quantify the 
amount of static order imposed by the confinement and 
have direct quantitative access to the influence of the 
set of frozen particles on the fluid structure to obtain 
bulk, equilibrium, many- body information not contained 
in pair correlations [7HlH Il7| . 

^From Fig. [2l we also conclude that measuring the 
evolution of the static overlap Qca is more difficult 
than previously thought [1, [lO, EH, because the dy- 
namics slows down considerably with increasing confine- 
ment [1^1 . Whereas the bulk dynamics at T = 8 corre- 
sponds to a moderately viscous state, there exists in all 
four geometries a maximal confinement above which Q(t) 
does not reach its long-time limit in the time window of 
our simulations, and hence large values of Qoo cannot be 
accessed. This difficulty gets even more pronounced at 
lower T [l7|, meaning that the measurement of point- 
to-set correlations closer to T c is, at present, extremely 
challenging. 

In Fig.[5J we present the evolution of the average over- 
lap Qoo(C) as a function of the confining length £ in all 
geometries (a), and the evolution of the corresponding 
relaxation times r(£) (b), defined from the time decay of 
Qself(^) to the value 1/e. In the latter figure we normalize 
the data by the bulk value Too = t(£ —> oo). The static 
profiles confirm that the average overlap becomes increas- 
ingly non-random by increasing the confinement since at 
a given confinement length we find that Qoo increases 
for the sequence pinned, wall, sandwich, box. Note that 
for the first three geometies the values of Qoo are simi- 
lar, whereas the ones for the box are significanly larger. 
This indicates that the effect of confinement is highly 
non-linear. In cases (b) and (d) we find that the large 
£ decay of Qoo is compatible with an exponential decay, 
Qoo(C) - Qrand ~ exp(-£/f stat ), which defines a PTS 
static lengthscale, £ s tat- For geometry (d) this depen- 
dence holds for temperatures as low as T — 5.0 < T c p7| . 
A compressed exponential decay was reported for a closed 
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FIG. 4: a) Dependence of the overlap Qoo on the confining 
length for all geometries at T = 8 (symbols) fitted with an 
exponential decay at long distances (straight lines), b) Re- 
laxation time of the self overlap normalized by its bulk value. 



cavity in Ref. [10|, which is unfortunately the geometry 
for which our data is the most limited. For randomly 
pinned particles, one expects the non-random part of the 
overlap to scale as c = 1/f 3 for large £, and an exponen- 
tial decay is not expected in that case, which requires a 
separate analysis (see also (201). 

Regarding dynamics, we find that decreasing £ at con- 
stant T leads to a strong slowing down of the dynam- 
ics for both self and collective quantities. In absolute 
values, t(£) strongly increases with the number of con- 
fining walls which shows that relaxation is non-linearly 
suppressed by adding more constraints. For a given con- 
fining length, pinned particles have of course much less 
impact on the dynamics as a single particle has obvi- 
ously less effect than an entire wall. For the cavity (a), 
the slowing down is in fact so dramatic that the range of 
£ where Qoo can be measured is very small: The overlap 
is too small when £ is large, but dynamics is much too 
slow when £ is small, which only leaves a narrow range 
to measure the static profile. We conclude that measur- 
ing PTS correlations in a closed cavity is in practice a 
difficult task in the interesting supercooled regime. 

For a single wall, r(£) can be followed down to £ — > 0, 
while for the other geometries small values of £ are not 
accessible due to a much stronger slowing down. All dy- 
namic profiles in Fig. HJd are well described by an expo- 
nential decay, ln(r/roo) ~ exp(— £/£dyn) which directly 



allows one to extract a dynamic correlation lengthscale 
£dyn [3 [HI • In Ref. [13] we have explored its tempera- 
ture dependence for a single wall and discussed how £dyn 
relates to previous measures of dynamic lengthscales. 
^From Fig. [4}d we see that the slowing down of the dy- 
namics is the least pronounced for the pinned geometry. 
Therefore we suggest that this geometry might be best 
suited best for the investigation of the T— dependence of 
static correlations. 

By comparing the two panels in Fig. 2J it is obvious 
that for all geometries the dynamics seems to be af- 
fected over a broader range of confinement than statics, 
which suggests that, generically, £dyn > £stat- In fact we 
have frequently found parameters for which clear devia- 
tions from bulk dynamics are observed while the static 
overlap is still random, suggesting a clear decoupling be- 
tween static and dynamic correlations. This decoupling 
is strong for a single wall and becomes less pronounced 
for randomly pinned particles. Our study thus seems to 
confirm [H, [lTj that static correlation length scales are 
generically decoupled from dynamic ones and smaller, at 
least over the range of temperature we can explore nu- 
merically. 

We finally discuss our results in a broader context. 
Theoretical progress on the glass problem is slowed by 
the lack of an obvious structural indicator to distinguish a 
glass from a fluid. It has only recently been realized that, 
in the framework of RFOT theory, the 'ideal glass' phase 
below the Kauzmann transition, Tr-, is characterized by 
an infinite static PTS lengthscale [1, Q . This implies for 
instance that below Tr-, freezing a semi- infinite space as 
in Fig. [1J1 determines the position of the particles in the 
entire z > space [6]. Spatially extended static profiles, 
such as shown in Fig. [4j thus uniquely characterize flu- 
ids approaching the glass transition [T3, E3 ; an d do not 
provide valuable information for liquids at high temper- 
atures. 

Within RFOT theory, a static length £rfot (T) 
emerges from a spatial interpretation of the (mean-field) 
concept of metastable states [l2j]. This 'mosaic' length 
scale plays a role similar to that of a nucleation length 
scale in first order transitions, because it is set by the 
competition between the entropic gain of exploring dif- 
ferent states and the energy cost of having interfaces be- 
tween them [l2j]. This should be reflected, in the closed 
cavity of Fig. QJi, as a crossover between a small Qoo 
when d > £rfot, to a large Qoo when d < £rfot, since 
for decreasing confinement the surface tension eventually 
dominates |8| . A qualitatively similar crossover holds for 
geometries (b, c) as well, but not for a single wall in (d) 
where bulk behaviour with — Q ra nd is recovered far 
from the wall and an interface must always be present. 

In contrast to closed cavities, the crossover set by 
£rfot involves in geometries (b) and (c) a number of 
particles that diverges in the thermodynamic limit. This 
'crossover' should therefore more properly be described 
as a genuine freezing transition towards an ideal glass 
phase where is large and density fluctuations do not 
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relax (see also [22[). A glass phase can then be ap- 
proached either by decreasing T in the bulk, or by in- 
creasing the confinement at constant T, which opens ex- 
citing perspectives to study the glass transition. We are 
presently pursuing the exploration of the glass phase ob- 
tained at large pinning density, but the results go much 
beyond the theme of the present article and will be pre- 
sented elsewhere. Our results also motivate further anal- 
ysis of the phase diagram of viscous liquids in confined 
geometries (b, c). Analytical calculations for randomly 
pinned particles exist for hard sphere systems within 
mode-coupling theory [23| , while thepresent results have 
motivated both an RFOT analysis 22[ and some numer- 



ical investigation [20(. Future work should characterize 
and compare in more detail the temperature evolution 
of the static and dynamic lengthscales introduced in this 
work beyond the case of the single wall studied in [TtJ • 
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